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Quantum mechanics in the vicinity of black holes is a fascinating field of theoretical physics. It 
involves both general relativity and particle physics, opening new eras to establish the principles of 
unified theories. In this article, we show that quantum bound states with no classical equivalent 
-as can easily be seen at the dominant monopolar order- should be formed around black holes 
for massive scalar particles. We qualitatively investigate some important physical consequences, in 
particular for the Hawking evaporation mechanism and the associated greybody factors. 
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' Black holes are extreme objects whose study is very rich and relics on different branches of physics. Among the 
most important and fundamental phenomena of black hole physics are the Hawking evaporation mechanism -see [l[ 
for a review- and the existence of quasi-normal modes (QNMs) -see for a review. On the one hand, the evaporation 
phenomenon reveals the profound links between gravity and thermodynamics. Its study is both extremely fruitful in 
CN| ■ itself and because of the quantum gravitational effects expected to occur during the last stages of the evaporation, 
when the semi-classical approach breaks down. On the other hand, quasi-normal modes are of particular importance 
because they revealed black holes to be stable under perturbations and represent a key ingredient in the computation 
' of gravitational wave signals. 

This article focuses on the investigation of bound states for massive particles emitted by black holes. The problem 
f^i ■ is mathematically quite close to the investigation of massive QNMs [3|, as in both cases the point is to find the 
\ characteristic complex frequencies allowing a massive (scalar) field to propagate in a black hole background, while 
satisfying given boundary conditions at the black hole horizon and at spatial infinity (see, e.g., [H). Nevertheless 
"p^ ' because those boundary conditions arc not the same, the physical meaning of bound states and quasi-normal modes 
+^ . are very different and correspond to different energy ranges. Although bound states are well known to exist in quite 
O ' a lot of classical and quantum systems, this work points out their specific existence around black holes, even at the 
D , monopolar order, and investigates some important consequences. In most of the literature, bound states refer both to 
' particles orbiting the black hole or back scattered to the black hole. In this work, we are interested in quanta trapped 
| between the finite black hole potential barrier and the infinitely thick well which prevents the particles from reaching 
infinity if their energy is low enough. As will be demonstrated hereafter, those orbiting states are not strictly stable 
but can be characterized by finite lifetimes. 
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In Section II, we first determine the conditions for such states to exist and compute their energy spectra and 
lifetimes at the WKB order. Thanks to a simple toy model, we estimate in Section III the average number of trapped 
particles as a function of their characteristic energies. Those trapped massive particles -whose greybody factors must 
be computed numerically- will also inevitably modify the Hawking radiation spectrum at infinity, which is computed 
in this Section. Some conclusions are finally given together with perspectives. 



II. ORBITING BOUND STATES 



The investigation of orbiting quantum bound states around Schwarzschild black holes requires one to solve rela- 
tivistic quantum mechanical equations in a curved background, while taking into account a non-vanishing mass. To 
show that those states do exist, the Klein-Gordon equation in a Schwarzschild background will be shown to exhibit a 
radial potential containing a local well for given ranges of black hole horizon radii and particle masses. 
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A. Conditions for such states 



The Klein-Gordon equation of motion for a scalar field $ with mass fi in a space-time with metric g a p can be 
expressed as 

^fy*) + p 2 $ = 0. (1) 

Writing <f> = e~ iut Y^ l (9, (p)R(r) to split the temporal, angular and radial parts of the field (where Y^ n arc the spherical 
harmonics), the radial function R(r) obeys, in a 4— dimensional Schwarzschild background, 



h(r) d , . , 2 d 2 , , . (£{£+l) 2 
_Y —h(r r 2 — + c^ 2 - fc(r) -^-^ + /i 2 
r^ dr dr \ r 



R(r) = 0, (2) 



where h(r) is defined by the metric ds 2 — h(r)dt 2 — dr 2 /h(r) — r 2 dVl 2 (see, e.g., [j| and references therein for a 
description of the general techniques associated with quantum fields in a Schwarzschild spacctime used throughout 
this paper). Under the change of variables r — > r* and i?(r) — > C/(r) where is the tortoise coordinate (such that 
dr* = dr/h(r)) and f/(r) = rR(r), this equation takes a Schrodinger-like form 

^ + (. 2 -V?(r))[/ = 0, (3) 



with a potential 



where rjj stands for the Schwarzschild radius and £ for the angular quantum number. The usual quantum mechanical 
techniques can therefore be employed in the tortoise coordinate system. The Chandrasekhar convention is used 
hereafter : the last term of equation ([3]) is interpreted as the squared potential so as to recover the standard Hamilton- 
Jacobi equation. On Fig. [H V 2 (r) is shown for three different values of the mass (/it = {0, V0.1, V0.4} D"^ 1 ]) and 
two values of the angular momentum {£ = and £ = 1). Depending on fi and £, it can be seen that a local minimum, 
potentially leading to a bound state, eventually appears. The existence of a potential well depends on the roots of 

the algebraic equation = : 

r H ^ 2 r 3 - 2£(£ + l)r 2 - 3r H (l - £(£ + l))r + Ar 2 H = 0. (5) 
Two roots above rjj exist if the mass \i is lower than a critical value n+{£) given by 
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with 



/4 = I - 27J i + \/ 729J i 2 + 432J 2 ) , (6) 



Ji = £ 3 (£+lf +£ 2 (£+l) 2 -?{£+l)-l, (7) 
J 2 = £ 2 (£ + l) 2 [9 + U£(£ + 1) + 9£ 2 (£ + l) 2 ] . (8) 

In the monopolar case (i.e. £ = 0), it takes the simple value n+ = \r~^ . This allows one to understand easily the 
general behavior of the potential. Let us consider a particle with a given mass [i\ and a black hole with a horizon 
radius r\. The mass of the particle may be smaller than l/2ri. In this case, the local barrier and the potential well 
do exist for all the values of the orbital quantum number, and particles can be potentially trapped with any angular 
momentum. On the other hand, if [i\ is greater than l/2ri, then the ^-domain has to be divided into two subclasses. 
A critical value of the orbital quantum number, denoted £\, can then be defined by /^+(£i) < fx < H+(£i + 1)- For all 
the partial waves with an orbital quantum number smaller or equal to this critical value, the potential appears as a 
monotonically increasing function of r from V^( r i) = to Vf(+oo) = However, all the partial waves with £ > £i 
will face the local barrier as well as the well potential. 

The detailed shape of the potential is also determined by another critical mass, hereafter called ^-(£), which defines 
the relative height of the potential barrier close to the horizon when compared with the mass of the particle. If the 
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FIG. 1: Square of the potential as a function of the radial coordinate for three values of the mass fi and two values of the 
angular quantum number £ (in each case, £ = for the lower curve and I = 1 for the upper curve). 



mass of the particle is higher than the barrier, there is no more a turning point for the system. The mass is the 
highest value of the potential if the equation V 2 (r) = /j, 2 has no root above ru- The cubic form 



satisfies this criterion if /i > /!_(£), which is then given by 

1 



%2 l-^+l) M r g _ Q 



27r% 



-Lj. + JLj + 27L 2 , 



with 



Li 
L 2 



= 2£ 3 {e + l) 3 + M 2 (£ + l) 2 - M{£ + 1) 
= ^+1) 2 [^(£+1) + 1] 2 . 



(9) 



(10) 



(11) 
(12) 



In the monopolar case, it becomes /i_ = ~^ t h ■ Table U provides the values of the critical masses in units of r^ 1 . 



t 




M+ [th 1 ] 




1 


•2 

0.794 


i 

2 

0.931 


2 


1.275 


1.480 


3 


1.768 


2.046 


4 


2.264 


2.617 


5 


2.761 


3.191 



TABLE I: Critical masses in units of r^ 1 for different angular quantum numbers. 

The potential is always zero at r = th and tends to [i for r — » oo. i) If (i < fj,— (£), the potential reaches a maximum 
higher than \i and then reaches a minimum; if ii) fi = /i_(£), the maximum is exactly equal to \i and a minimum 
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also appears; if iii) H-(Jt) < fJ, < the potential reaches a maximum lower than [i and still admits a minimum, 

whereas if iv) fi > fJ.+ (£) the potential is a monotonically increasing function of r. This behavior is illustrated for 
£ = on Fig. [2] for masses corresponding to those four specific cases. As a direct consequence of the wave description 
of particles in quantum mechanics, it can be seen that bound states -due to the local minimum- can appear at the 
monopolar order, with no classical equivalent. Although particles can of course be classically trapped around a black 
hole, no state without angular momentum can be found if the quantum behavior is not taken into account. Just as 
in classical mechanics, it is the non- vanishing mass that leads to the trapping (even at the quantum level), but with 
specific quantum features, lying in the quantization of the angular momentum, and furthermore allowing for trapping 
even without any angular momentum. Furthermore, whatever the mass of the particle, a bound state will appear for 
high enough multipolar orders so that fi < fJ,+ (£). This makes this phenomenon of "particle trapping" quite generic. 
Those bound states are described by quasi-stationary quantum states that cannot reach spatial infinity but can still 
make a transition to the black hole by tunneling through the gravitational barrier. 




FIG. 2: Squared monopolar potential for five particle masses : i) /i 2 = 0.1 r H 2 corresponding to fj, < /u_(0), ii) /i = ^_(0), iii) 
/i 2 = 0.2 r^ 2 corresponding to /U-(O) < /x < /U+(0), iv) fi = /i+(0) and v) /j 2 — 0.4 r^ 2 corresponding to fi > n+(0). 

As the qualitative features can easily be understood at the monopolar order, this particular value of the quantum 
angular momentum is now assumed. When \i < the positions of the potential barrier (r_) and of the local 
minimum (r + ) can be analytically determined to be 



1 

r± = - 
/i 



cos (6»)±v / 3sin(6») 



(13) 



with 



9 = — arctan 
3 



^1 -1 



(14) 



The asymptotic behavior is in agreement with the monopolar potential for a massless particle : 



lim r. 



+oo , lim r_ = — rj{ 
fj.^o 3 



lim r± = 2th- 



(15) 



The latter case, \x — > corresponds to a saddle point at r = 2r^, which represents the degeneracy of the maximum 
and minimum of the potential. The most important properties of the monopolar potential are summarized in the 
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various graphs of Figure [H As long as the particle's mass is smaller than the potential displays a local well and 
this feature may lead to the existence of quasi-stationary states localized in the well. 




Particle's mass [r '] Particle's mass [r '] 



FIG. 3: Upper left : position of the local minimum of the potential Vmin, in units of th , as a function of the mass of the particle. 
Upper right : ratio of the local minimum over the mass of the particle, Vmm//i, as a function of the mass of the particle. Lower 
left : position of the local maximum of the potential V max , m units of rH, as a function of the mass of the particle. Upper 
right : ratio of the local minimum over the local maximum, Knm/Knax, as a function of the mass of the particle. 



B. Complex frequency spectrum : a WKB analysis 

Even when there is a well in the potential, this does not guarantee orbiting bound states to exist. To ensure 
the stability of those states, the local well has to be deep enough, because of the zero-point energy associated with 
quantum systems. Under the reasonable assumption that the shape of the potential is well approximated by a second 
order expression, the dynamics is similar to an harmonic oscillator with a frequency u> given by the curvature of the 
potential around its minimum : uj = ^ d 2 V 2 / dr- k 2 \ r+ . This curvature must be evaluated as a function of the tortoise 

coordinate, because this is the coordinate system in which the radial part of the Klein-Gordon equation is of the 
Schrodinger type. The zero-point energy is simply given by wo = The ratios of this approximated computation 
to the mass of the particle and to the maximum of the potential are plotted on Fig. [4] Each time a minimum does 
appear, the zero-point energy remains smaller than the mass and the gravitational potential barrier, allowing bound 
states to exist. 

Because the left part of the potential well is not an infinite barrier, the orbiting bound states are describes by quasi- 
stationary states associated with resonances in the density of states. Those resonances are characterized by their 
complex frequencies, the real part corresponding to the energy position and the imaginary part to the bandwidth. 
The spectrum of complex frequencies can be determined at the WKB order using the techniques developed in 
The spectrum will be infinite if fi < /i_ and finite if /i_ < fx < The Bohr-Sommerfeld rule, whose validity in 

a relativistic framework was established in [9(, reads V n € N, / ^Juj 2 -V 2 {r) — — = \ n-\ — 7r, the highest 

frequency allowed for a bound state being [i for [i < /x_ and V(r_) for /i > In the latter case, the left-hand side 
integral is clearly finished and there exist n max states. If fj. < [i- , the upper bound of the integral is infinite when 
u) = fi and, the function to be integrated being proportional to r -1 ' 2 near +oo, the spectrum is expected to be infinite. 
With the appropriate change of variables in the cubic equation giving the turning points, it can be shown that the 
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FIG. 4: Ratio of the zero-point energy ujq and the mass jj, and ratio of ujq and the maximum of the potential V m ax, both as a 
function of the mass of the particle in units of . 

normal frequency spectrum of the resonances depends only on [ir^ and can be written as {uj n } ne f} = \ ^^"^ \ 

Table UTl gives some normal frequencies and the associated bandwidths, as numerically obtained at the semi-classical 
order, following (Io| to evaluate the tunnel probability. The quantity N eq in this table will be explained in the next 
section. 



C. Qualitative features of the bound states and comparison with quasi- normal modes 



Before investigating the consequences of bound states for Hawking radiation, some qualitative features of those 
states, related to QNMs, and the existence of a halo of trapped particles around black holes, are briefly discussed in 
this section. 

As underlined in the computations of bound states and QNMs are very similar. The latter are given by the 
pure outgoing boundary conditions [f| while the former require evanescent waves at spatial infinity and ingoing modes 
at the black hole event horizon. However, there are no other fundamental physical links between QNMs and bound 
states. First of all, quasi-normal modes are related with the energy carried out by perturbations of black holes whereas 
bound states correspond to quantum states that cannot escape from the gravitational attraction. In addition^, the 
mass fj, of the field has to be smaller than H-(£) for QNMs to exist at the £th multipolar order. Moreover, because 
quasi-normal modes correspond to resonances near the peak of the potential, their energies have to be greater than 
the mass of the field. Those characteristics are direct consequences of the pure outgoing boundary conditions. The 
situation is fundamentally different for orbiting bound states. The condition on fi is less restrictive : \i has to be 
smaller than /j,+ (£) > /J.~(£) for orbiting bound states to exist at the £th multipolar order. Furthermore, the energy 
of those states is lower than the mass of the particle, due to the gravitational binding energy. 

In addition, as long as /i < /i+(0), bound states exist for £ = and a spherical halo of quanta "orbiting" the black 
hole can be expected. If the mass is between fi+(£ — 1) and bound states will exhibit an angular distribution 

dominated by the lowest multipolar order allowing for a minimum in the potential, that is with a distribution roughly 
given by (6, ip). Clearly, the mass of the field has to be close to the critical masses /z± for the orbiting bound states 
to have a substantial influence. In particular, if fi <C ^ + (0), the local well is so tiny at any multipolar order that 
particles can be treated as massless. So as to fix the orders of magnitude, table UTTl gives the masses and temperatures 
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= 0.25 x n+ 




n 





1 


2 


3 


4 




1.24.KT 1 


1.25.10 -1 


1.25.10 -1 


1.25.10 -1 


1.25.10- 1 


Tr H 


1.49. 10~ 2 


1.49. 10" 2 


1.50.10~ 2 


1.50.10~ 2 


1.50.10" 2 


N eq r H l 


2.15 


2.09 


2.09 


2.09 


2.09 





H = fx- 


n 





1 


2 


3 


4 


LOTH 


3.75.10 -1 


3.83.10 -1 


3.84.10 -1 


3.84.10" 1 


3.85.10 -1 




1.68.10- 1 


1.87.10- 1 


1.90.10- 1 


1.91.10 -1 


1.92.10" 1 


N eq r H l 


0.024 


0.021 


0.021 


0.021 


0.020 






fj, = 0.8 x /u+ 


n 





1 


2 


3 


4 


UJT H 


3.88.10- 1 


/ 


/ 


/ 


/ 


Tr H 


1.91.10 -1 


/ 


/ 


/ 


/ 


N eq r H L 


0.019 


/ 


/ 


/ 


/ 



TABLE II: First and second rows below the integer n : Spectrum of normal frequencies uj and bandwidths Y evaluated at the 
WKB order in units of r]j . Third row : Number of trapped particles N eq in the different energy levels once the equilibrium 
regime is reached (see the text). Three different masses are considered : /j, < /j,- , [i — [i- et [i > /i-. 

black holes should have so that the effects studied in this article become important for some standard model particles. 
Although they are not spinless -therefore requiring one to investigate the master equation for fermions and gauge 
bosons- the main qualitative features can safely be inferred from the scalar case. 



particle 
(mass) 


M = M- 
Mbh [kg] 


for I = 
T H [MeV] 


fi = n+ f 
M BH [kg] 


or I = 
T H [MeV] 


electron 
(511 keV) 


9 x 10 la 


0.1 


1.2 x 10 14 


0.08 


muon 
(105 MeV) 


5 x 10 11 


21.7 


6 x 10 11 


16.7 


charm 

(~ 1.3 GeV) 


4 x 10 iu 


248 


5 x 10 iu 


190 


top 

(171 GeV) 


2.5 x 10 s 


36 x 10 d 


3.7 x 10 s 


28 x 10 d 


W boson 
(80 GeV) 


6 x 10 8 


16 x 10 a 


8 x 10 8 


12 x 10 a 



TABLE III: Masses and temperatures expected for a black hole to have the critical masses fi± close to the masses of some 
standard model particles. 

Most very massive particles (like the top quark or W boson) are unstable and will decay before reaching infinity 
or creating a substantial halo around the hole. This point should, however, not prevent us from considering bound 
states associated with lighter stable particles, as the relevant parameter for the bound states to exist is not the mass 
of the field in itself, but the hierarchy between this mass and the mass of the hole. As far as the field is not strictly 
masslcss, there will be a black hole masse range where ^ becomes close to fx+ and fj,-, making the trapping effective. 

An important issue to address is related to the time stability of those bound states. By integrating the Hawking 
instantaneous energy spectrum and summing over all the degrees of freedom of the standard model of particle physics, 
it is easy to show that the mean time between the emission of two quanta is, in Planck units, At « 100M. If we 
consider, e.g., electrons emitted by a black hole such that /i « this leads to At « 10 _20 s; this remains tiny when 
compared to the time-scale of the black hole evolution. In other words, as the energy carried out from the hole by 
each emitted particle (E w T = l/(8nM)) is much smaller than the mass of the hole, many quanta can be emitted 
and trapped before the structure of the potential will have substantially changed. This is, of course, not true anymore 
when T ~ M, that is to say in the Planck regime in which the semi-classical approximation breaks down anyway. 
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III. CONSEQUENCES FOR THE HAWKING RADIATION 



A. Trapped particles 

The existence of bound states will play an important role in the Hawking evaporation mechanism which was initially 
described in Q without taking into account this phenomenon. 

First, some low-energy particles will be trapped by the potential well. They will not reach infinity and the spectrum 
will be modified in a way quite similar to what could happen due to a QCD halo [!]• When the black hole evaporates, 
new bound states will appear each time the temperature becomes of the same order than the critical masses /i± 
associated with existing particles. At the first level of approximation, the number of particles within the well per 
unit of time and energy, dN/dt, at a given energy is given by Fi„ — F out with Fi n = uj /\r(e u t T ± 1)] representing the 
ingoing one due to the Hawking evaporation, and F out = N/t the outgoing one. The energy-dependent time constant 
t = 1/r simply corresponds to the lifetime of the bound state. This simple model for the number of particles trapped 
in a local well is realistic only in the case of heavy black holes which can be considered as stationary (the small amount 
of energy carried out by emitted particles being much smaller that the mass of the black hole). With Nq the initial 
amount of trapped particles, the time evolution is simply given by 

N( U , t) = N (.)e-^ + (l - e~^) . (16) 

which can be clearly understood as competition between particles "leaking" from the well to the hole and particles 
"filling" the well because of the Hawking radiation. Quickly the number of trapped particles reaches an equilibrium 
regime with N eq ~ ^^_^ / 1 Tg — j-j. It is worth noticing that once the equilibrium regime is reached, the number of 

trapped particles does not depend on r but only on the energy and on the temperature of the black hole. This 
behavior is due to the fact that the probability to cross the potential barrier is the same for incoming and outgoing 
particles. As a consequence, the equilibrium regime, given by the ratio of the number of particles scattering the off 
barrier from the left and the number of particles scattering off the barrier from the right, does not depend on the 
lifetime of the bound states. 

It is important to compare the characteristic time to reach this regime with the characteristic evolution time of the 
black hole. For monopolar bound states, which are dominant, and for particles with masses below the lifetime 
ranges between 10M and 100M. This time-scale is of the same order as the typical time-scale between two successive 
emissions of a particle , which, as already mentioned in this article, is much smaller than the evolution time-scale for the 
black hole itself. This makes meaningful our hypothesis of non-evolutionary black holes : as long as heavy black holes 
are considered, enough particles are emitted to reach the equilibrium regime without any substantial modifications of 
the properties of the black hole. The expected number of trapped particles when the equilibrium regime is reached is 
also given in Table HT1 in units of r^} 1 , per unit of energy and time. As can be seen from this table, the mean number 
of particles decreases as the mass of the particle increases, since higher energy bound states are involved. 



B. Greybody factors and radiation spectra 



The mass of the particle will also drastically modify the greybody factors that account for the non-trivial part 
(gravitational barrier and centrifugal potential) of the couplings between quantum fields and evaporating black holes. 

I ) have recently been computed 
I , rotating black holes [IH, , 
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m followed by [13|,11J : 
|7de-Sitter spacetime 



The greybody factors (whose detailed study began with 
in quite a lot of interesting situations : extradimensions [1| 

Gauss-Bonnet gravity [2(| [HJ etc., but up to now the masses of the emitted particles have mostly been ignored 
(although some good estimates were obtained in [HI, [22j). Figure O displays the absorption cross section numerically 
computed by solving the Klein-Gordon equation to evaluate the ingoing and outgoing amplitudes of the wave function 
at the horizon and at spatial infinity (see, e.g., [l7j for a detailed description of the method we have developed). It 
should be pointed out that those cross sections have been computed for particles with an energy greater than or equal 
to the mass. Particles with an energy smaller than the mass will face an infinitely thick potential barrier, precisely 
due to the mass term in the potential, and will never reach spatial infinity. It should also be noticed that, when 
the wavelength of the particle becomes infinite, the cross section diverges, potentially leading to an experimentally 
relevant enhancement of soft quanta. It can indeed be expressed as 



a g( UJ ) = 



ir{2£+ 1) 



1-4, 
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where k is the momentum (so that uj 2 = k 2 + /i 2 at spatial infinity) and \ Ag\ is the transmission coefficient. However, 
the transmission coefficient remains always smaller than or equal to unity and the divergence lies in geometrical issues. 
When w — > fx, the wave number of the particle tends to zero, whereas the transmission coefficient is non-vanishing, 
just because the potential barrier is not infinitely thick (the barrier precisely stops where the potential well starts 
and the potential tends asymptotically to /i with negative values). This leads to a k~ 2 divergence. The situation is 
clearly different for massless particles. In this case, the potential barrier tends asymptotically to zero with positive 
values, leading to \Ai\ 2 oc k 2 when u> — > 0. This feature prevents the IR divergence from occurring for massless 
particles. Those results are in agreement with analytical investigations in the IR regime, a g (k — > 0) ~ Air 2 (fir h) 3 / / k 2 , 
as obtained following [22| . 

In Fig[5J the flux at infinity emitted by a black hole is plotted when the masses of the emitted quanta are taken into 
account. As can be seen, this substantially modifies the usual picture both because of the intrinsic cutoff imposed 
by the mass and because of more subtle effects included in this analysis, like the selection induced on the allowed 
quantum multipolar orders of the outgoing particle. From Fig[6j it can be seen that the radiation flux decreases for 
higher masses. However, in the IR limit, this tendency should change, as can be seen using the analytical limit 



d 2 N 
dtdk 



(*-0) 



ATT 2 (pr H ) 3 
e^l T - 1 



obtained from [221 ] . The flux is non- vanishing only if the mass of the particle is non- vanishing, enhancing the emission 
of ultra-soft quanta, just as for Schwarzschild-de Sitter black holes [17j . 
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FIG. 5: Absorption cross section for massive scalar particles in units of nr H as a function of the energy measured at infinity. 



IV. CONCLUSION 



This study establishes the existence of new bound states around black holes which, at least at the dominant 
monopolar order, have no classical equivalent. Although such states are known to exist in other physical systems, 
this opens new perspectives to investigate the detailed features of the Hawking spectrum (with possible cosmological 
consequences related, e.g., to the primordial power spectrum -sec [23] for recent limits and |!24] for a review), the 
intricate shape of the grcybody factors and the propagation of massive quantum fields in the vicinity of a black hole. 
Not only could the phenomenology be revised as the spectra should be quantitatively modified but fruitful thought 
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experiments associated with light black holes should also take into account those states. The intricate problem of 
backreaction should however be addressed and requires an exhaustive study in itself. 




FIG. 6: Flux at infinity emitted by a black hole when taking into account the mass of the emitted particle. The different curves 
share the convention of Fig. [6] : fj, = [r^ ], y/0.1 [rjj 1 ], /J,- , \/0.2 [rjj 1 ], fi+ [r^j 1 ], V0.4 [rjj 1 ] from top to bottom. Left panel : 
Hawking's radiation spectrum, d 2 N/dtdk, as a function of the particle's momentum k. Right panel : Hawking's radiation 
spectrum, d 2 N/dtduj, as a function of the energy lo of the particle. 
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